Design and optimization of wall-climbing robot impeller by genetic algorithm based on computational fluid dynamics and kriging model

In recent years, wall-climbing robots have begun to replace manual work at heights to reduce economic losses and casualties caused by working at heights. This paper designs a negative pressure adsorption type wall-climbing robot and analyzes the internal fluid movement state of its negative pressure device and the force analysis of the robot when it is adsorbed and balanced. Furthermore, through the experimental prototype, the influence of wall material, robot pose, negative pressure cavity shape and sealing method on the adsorption performance of the wall-climbing robot is explored. The computational fluid dynamics simulation (CFD) simulation method and experimental results are used to verify each other, which proves the correctness of the simulation results. Based on the Kriging surrogate model, the functional relationship between the impeller blade outlet angle, the impeller inlet diameter, the number of blades as the design variables, the negative pressure as the dependent variable was established, and the genetic algorithm (GA) was used to optimize it. Compared with the original design, the optimized design results of impeller parameters have increased the negative pressure value from 3534.75 to 4491.19 Pa, an increase of 27.06%.

ability of wall-climbing robots, there are more and more researches on rotating-flow negative pressure adsorption devices. Zhao measured the suction characteristics of the suction unit and studied the influence of the wall roughness and shape on the adsorption performance 23 . Amakawa developed a wall-climbing robot attached to the fuselage of an aircraft based on negative pressure technology. Through basic experiments and calculation of the adsorption force, the relationship between the adsorption force generated by the centrifugal fan and the size of the adsorption part, the number of revolutions of the centrifugal fan, the leakage area, and the characteristics of the centrifugal fan and the influence of the curved surface on the suction force are revealed 24,25 . Computational fluid dynamics simulation (CFD) and other methods can be used to analyze, improve and optimize the adsorption device more conveniently 26 . Chen studied the blade vortex's theoretical model, modeled the rotating blade's height, analyzed the secondary flow of the blade by the CFD method, and proposed a blade design method based on the characteristics of the motor 27 . Fallah performed numerical calculations and CFD analysis on the wall-climbing robot designed by his vortex current adsorption to obtain the optimal data and designed a new negative pressure cavity structure 28 . The vortex flow device can form negative air pressure and is also suitable for other fluids. Zhu introduced an electrically-driven underwater whirlpool suction device. The experiments and CFD calculations show that the suction increases quadratically with the increase of speed. Moreover, as the gap increases, the suction first increases and then decreases 29 . Koo gave the static and aerodynamic models of the adsorption mechanism and analyzed them, which provided a basis for design and control 30 . Using a gas pressure sensor to monitor the pressure value in the negative pressure cavity and keeping the pressure in the negative pressure cavity constant can improve the stability of the wall-climbing robot when it moves 31 .
The main contributions of this paper are as follows: 1. In this paper, an experimental prototype of negative pressure adsorption wall-climbing robot is designed. The parameters that affect the adsorption performance of wall-climbing robots are studied, such as wall materials, robot pose, negative pressure cavity shape and sealing method. 2. A new optimization method is proposed for the impeller. The sample points required for the kriging model were obtained by CFD simulation. A surrogate model between impeller parameters and negative pressure values was established using the kriging model. The impeller parameters are optimized based on the genetic algorithm, which improves the adsorption capacity of the robot.
The rest of this article is organized as follows. In "Analysis of fluid state and adsorption conditions" section, the motion state of the fluid in the adsorption device of the wall-climbing robot is analyzed, and the adsorption condition equation is established. In "Adsorption method of wall-climbing robot" section, a prototype of a wall-climbing robot with negative pressure adsorption is designed. In "Prototype experiment and CFD simulation verification" section, the parameters that affect the adsorption performance of the wall-climbing robot are experimentally studied. The CFD method is used for simulation and mutual verification with experiments. In addition, in "Optimized design of impeller" section, the critical components of the impeller are optimized based on the kriging method and GA, and the optimization effect is verified by CFD simulation. Finally, the conclusion is given in "Conclusion" section.

Analysis of fluid state and adsorption conditions
Negative pressure adsorption means that the pressure in a part of the area is lower than the atmospheric pressure. The relatively high-pressure atmosphere flows to the negative pressure area, thereby generating a specific force. If the negative pressure is too small, it will not absorb, resulting in slippage or even falling; when the negative pressure is too enormous, it will make the robot challenging to move and increase the robot's power consumption. Figure 1 shows the components of the negative pressure mechanism of the wall-climbing robot. The motor drives the impeller to rotate at high speed to generate centrifugal force. The impeller draws the gas in the negative pressure cavity, thereby forming a negative pressure in the cavity. The function of the sealing device is to reduce air leakage.
Analysis of fluid state in negative pressure cavity. Analysis of the state of fluid movement in the negative pressure cavity. When analyzing the fluid between the sealing method and the wall gap, the Reynolds number Re must first be used to determine whether the fluid flow is laminar or turbulent. When Re <2320 , the fluid motion was laminar; then Re >2320 , the fluid motion was turbulent 32 .
The calculation formula of Reynolds number Re is: where v is the average flow velocity of the fluid; d k is the hydraulic diameter of the flow section; η is the kinematic viscosity of the fluid, and η=µ/ρ , µ is the dynamic viscosity of the fluid, ρ is the density of the fluid, at room temperature, η=1.57 × 10 −5 m 2 /s. According to the similar principle, the hydraulic diameter of non-standard structure pipeline can be calculated by the following formula 33 Among them, A is the cross-sectional area of the flow; S is the sum of the inner and outer perimeters of the flow. www.nature.com/scientificreports/ The negative pressure cavity structure of the negative pressure adsorption type wall-climbing robot has a rectangle and a circle, and the structures are shown in Fig. 2a and b, respectively.
As shown in Fig. 2a, for the rectangular negative pressure cavity, the hydraulic diameter d ks between the sealing device and the wall gap δ is: Among them, l 1 is the dimension of the negative pressure cavity in the y direction; l 2 is the dimension of the negative pressure cavity in the x direction; a is the side length of the sealing method.
Take a = 0.1l 1 = 0.1l 2 , then formula (3) can be simplified to d ks ≈ 0.2l 1 , considering that the wall-climbing robot is generally small in size, taking l 1 ≤ 0.5 m , then the rectangular negative pressure cavity d ks ≤ 0.1 m −1 , the flow velocity between the sealing method and the wall gap is very low, taking v ≤ 0.1 m/s . From Eq. (1), the Reynolds number Re s ≤ 637 < 2320 of the rectangular negative pressure cavity can be obtained. Therefore, the  www.nature.com/scientificreports/ flow of fluid between the sealing method of the rectangular negative pressure cavity and the gap in the wall is laminar flow.
In the same way, the hydraulic diameter d kc of the circular negative pressure cavity is: Therefore Re c ≤ 637 < 2320 , the flow between the sealing method of the fluid-shaped negative pressure cavity and the gap in the wall is also laminar.
Force analysis of fluid movement in negative pressure cavity. Since the air velocity entering the negative pressure cavity through the gap between the sealing method and the wall is shallow and laminar, the Navier-Stokes equation (N-S equation) can be obtained: Among them, ρ air is the density of air, at room temperature (20 °C), its value is 1.20 kg/m 3 ; u is the fluid velocity, and its components in the x , y and z directions are u x , u y and u z respectively, ∇ is the gradient symbol; p is the pressure between the sealing device and the wall gap δ ; µ air is the dynamic viscosity of the air, at room temperature (20 °C), its value is 0.0183 × 10 −3 Pa·s; f is the unit mass force, and its components in the x , y and z directions are f x , f y and f z respectively.
Assuming that all the air particles in the gap are moving in parallel and straight lines, the speed u has nothing to do with x , y direction, that is u = u(z) 34 .
Choose the y axis direction as the air movement direction, only u y non-zero and ∂u/∂y = 0 , u x and u z are both 0, the N-S equation in the and directions is Then the pressure p has nothing to do with x , z , that is p = p y ; Assuming that the air in the gap is not compressible, there are: Substituting Eqs. (6) and (7) into Eq. (5), we can get: Performing quadratic integration on formula (8) can be obtained: Among them, C and D are integral constants, and the boundary conditions are: Substituting Eq. (10) into Eq. (9), we can get: The air flow q into the negative pressure cavity through the gap is Among them, p 0 is the standard atmospheric pressure value, which is 101.325 kPa; p 1 is the internal pressure of the negative pressure cavity of the wall-climbing robot, and the pressure distribution in the gap is shown in Fig. 3.
It can be seen from Fig. 3 that the air pressure outside the negative pressure cavity is p 0 , and the pressure becomes p 1 , after passing through the gap, resulting in a linear pressure distribution with a pressure difference of p = p 0 − p 1 , then the resultant force F δ of the atmospheric pressure acting on the sealing method is: µ air dp dy (9) u = 1 µ air dp dy The result is: Adsorption conditions of wall-climbing robot. It is assumed that the supporting force of the wall acting on the robot is borne by the four rollers on both sides of the track, yoz is the wall, the z-axis is the forward direction of the robot, and the x-axis is perpendicular to the wall. Figure 4 shows the force analysis diagram of the wall-climbing robot in any posture. Among them, N 1 , N 2 , N 3 and N 4 are the pressure of the wall on the left front wheel, left rear wheel, right front wheel and right rear wheel of the wall-climbing robot respectively; µ 1 is the static friction coefficient between the wall and each wheel of the wall-climbing robot; G is the gravity of the wall-climbing robot, and its magnitude is G = mg ; m is the mass of the wall-climbing robot; g is the acceleration of gravity; H is the distance between the center of gravity of the wall-climbing robot and the wall; L is the distance between the front and rear wheels of the robot; B is the gauge of the two crawlers of the wallclimbing robot; θ is the angle between the gravity and the orientation of the body.
The robot is balanced by forces in the x and z directions and moments around the two outermost rollers. The dynamic equation of the wall-climbing robot in any posture is: From Eq. (17), it can be seen that under the condition of any pose of the wall-climbing robot body, the differential pressure F p is: In order to ensure the normal work of the wall-climbing robot, it needs N 3 to be greater than 0, which can be obtained from Eq. (18). Under the condition of any posture of the robot body, the adsorption condition of the wall-climbing robot is:

Adsorption method of wall-climbing robot
Design of wall-climbing robot and its impeller model. Table 1 shows the basic parameters of the designed wall-climbing robot.
We need to design the appropriate blade size according to the load requirements of the robot. Therefore, the load quality of the robot in Table 1 refer to the pre-assumed impeller design conditions. The structure of the impeller plays a vital role in the size of the negative pressure. The main parameters of the impeller include impeller inlet diameter D 0 , blade inlet diameter D 1 , impeller outlet diameter D 2 , blade inlet size b 1 , blade outlet size b 2 , blade inlet geometric angle β 1A and blade outlet geometric angle β 2A 35 , as shown in Fig. 5. According to the minimum adsorption pressure under the adsorption equilibrium condition and the parameter design method of the impeller, the main parameters of the impeller can be obtained as shown in Table 2 36 , the three-dimensional model of the impeller established according to the impeller parameters is shown in Fig. 6a, and the wall-climbing robot model is shown in Fig. 6b. The closed impeller consists of a wheel cover, blades, and a wheel disc. The closed structure reduces internal leakage losses, so the working efficiency is higher than that of open and semi-open impellers, suitable for high-pressure situations.
Experimental prototype of wall-climbing robot. The actual object designed according to the abovementioned three-dimensional model parameters of the wall-climbing robot is shown in Fig. 7. The support www.nature.com/scientificreports/ structure between the impeller and the motor is made of high-strength nylon glass fibre 3D printing, which can support the motor and prevent the impeller from eccentricity during operation. The ground clearance of the negative pressure cavity can be adjusted by replacing the studs of different lengths in the crawler structure.
The influence of different ground clearance of the negative pressure cavity on the adsorption performance of the wall-climbing robot can be adjusted be analyzed. The main structure of the negative pressure cavity is made of high-strength acrylic plate material, which ensures the structural strength of the negative pressure cavity and facilitates manufacturing. The sealing device uses the sealing top produced by the 3 M company, which can reduce air leakage and achieve stable air pressure in the negative pressure cavity. The air pressure sensor is arranged for real-time measurement of the pressure inside the negative pressure cavity.

Prototype experiment and CFD simulation verification
Prototype experiment and analysis of its results. In order to control variables and avoid accidental errors, the speed of the driving motor of the impeller during the experiment is 8000 r/min. Furthermore, all the working condition experiments are repeated three times, and the final experimental value obtained is the average of the three experimental results.
The effect of wall material. In order to explore the adsorption performance of the wall-climbing robot when adsorbing on different wall materials, three different wall materials were used to conduct adsorption experiments. As shown in Fig. 8, the robot is adsorbed on the lime wall, the PVC (Poly Vinyl Chloride) wall and the ordinary corrugated box wall. The walls of lime, PVC and corrugated boxes have obvious grades in terms of hardness and surface roughness, and can have obvious contrast effects. In addition, it is easier to experiment with walls of these three materials in a laboratory environment. The experiment measured that the negative pressure value on the lime wall is 3351.24 Pa, the negative pressure value on the PVC wall surface is 3880.0 Pa, and the negative pressure value on the ordinary corrugated box wall surface is 4404.36 Pa. When the suction device of the robot is working, under the action of negative pressure, both the PVC material wall and the ordinary corrugated box wall are deformed. The deformation of the corrugated paper wall is more significant than that of the PVC material surface, which causes the gap between the negative pressure cavity and the wall to shrink. The air leakage between the gaps is reduced, thereby increasing the negative pressure value. At the same time, the lime wall surface has almost no deformation, and the gap between the negative pressure cavity and the wall surface has no noticeable change, so the negative pressure value is the smallest. When the wall deformation is not considered, the wall material has little effect on the adsorption performance.
The influence of robot orientations. Figure 9 shows the experimental results when the robot is adsorbed on the same lime wall in three different orientations, including three different robot orientations: vertical, inclined 45°, and horizontal. According to experiments, the negative pressure values of the wall-climbing robot in the three states are 3531.24 Pa, 3557.76 Pa, and 3578.16 Pa, respectively. The difference between the three experimental conditions is negligible. It can be considered that the posture of the wall-climbing robot adsorbed on the wall is opposite. The negative pressure value has no effect.  The influence of the shape of the negative pressure cavity. To compare the effects of different negative pressure cavity shapes on the adsorption performance, a circular negative pressure cavity (0.05 m 2 ) with the same area as the rectangular negative pressure cavity was used for adsorption performance experiments. Figure 10 shows the experiment of a wall-climbing robot in which the circular negative pressure cavity is adsorbed on the lime wall.
The experiments were carried out with three robot orientations consistent with the rectangular negative pressure device to control the variables to reduce the influence of other factors on adsorption. Table 3 lists the negative pressure values of two different cavity types. It can be seen that compared with the rectangular negative pressure cavity, the performance of the circular negative pressure cavity is better because the gas flowability of the circular negative pressure cavity is better.
Influence of sealing device. Two different negative pressure cavity sealing methods are shown in Fig. 11. The sealing method is shown in Fig. 11a used a sealing skirt to surround the bottom of the negative pressure cavity. The sealing method is shown in Fig. 11b is based on the first seal and then adds a layer to the bottom of the negative pressure cavity. The sealing skirt reduces the ground clearance of the robot. Two sealing methods are used to conduct experiments with the robot body upright. The experimental results are: the negative pressure value of the first sealing method is 3531.24 Pa, and the negative pressure value of the second sealing method is 3720.96 Pa. It can be seen that the second sealing method reduces the air leakage flow when the robot is adsorb-   www.nature.com/scientificreports/ ing and increases the negative pressure value by 5.4%. Therefore, a reasonable sealing structure is a key to the regular operation of the wall-climbing robot's adsorption device.
In summary, when designing the suction device of the wall-climbing robot, it is essential to consider the shape of the negative pressure cavity and the structure of the sealing method to improve the suction performance of the wall-climbing robot.

CFD simulation and verification. CFD setup and simulation.
In order to compare and analyze the experimental results, the CFD method is used to simulate and analyze the fluid movement inside the negative pressure device when the wall-climbing robot is working. The grid division of the fluid domain of the negative pressure device is shown in Fig. 12. The fluid domain mesh of the negative pressure device is mainly composed of tetrahedral and hexahedral mesh. The tetrahedral mesh is primarily located in the impeller fluid domain mesh. The hexahedral mesh is mainly located in the fluid domain mesh in the negative pressure cavity. Figure 13 shows the simulated stress cloud diagram of the impeller and the negative pressure cavity. Figure 13a and b, respectively, are the static pressure cloud diagram and the absolute pressure cloud diagram of the negative pressure cavity. The two stress distributions are the same because the static pressure p sta in Fluent refers to the difference between the absolute pressure p abs and the operating pressure p ope . The operating pressure set in the simulation is the standard atmospheric pressure 101,325 Pa. The pressure value in the static pressure cloud chart is negative, as shown in Fig. 13a. The negative pressure value is about 3500 Pa. The static pressure cloud diagram of the impeller fluid domain is shown in Fig. 13c. The maximum negative pressure is 19.8 Pa, near the outlet of the impeller fluid domain; the minimum is − 4183.3 Pa, near the impeller inlet. The dynamic pressure   Fig. 13d. The dynamic pressure p dyn is a physical quantity related to the speed, and its specific expression is:  20), it can be seen that the dynamic pressure is proportional to the quadratic of the fluid velocity. Because the fluid velocity in the negative pressure cavity is relatively slow, and the overall dynamic pressure is relatively small.
The impeller's dynamic pressure cloud diagram and velocity cloud diagram are shown in Fig. 13e and f, respectively. They are mainly divided into the impeller inlet area, the impeller outlet area, and the transition area between the two. The inlet area is characterized by low dynamic pressure and low flow velocity. In the simulation, the inlet dynamic pressure and velocity are almost 0; the outlet area has high dynamic pressure and high flow velocity. The maximum outlet dynamic pressure is 3380 Pa, and the outlet velocity is 83.1 m/s. As shown in Fig. 13, the transition area of the impeller is between the inlet and the outlet. The dynamic pressure and velocity values in the transition region are not constant. The dynamic pressure value in the transition zone is between the inlet and outlet. The same goes for the velocity values in the transition area.
The total pressure cloud diagram of the fluid area of the negative pressure device is shown in Fig. 13g and h, respectively. The calculation formula of the total pressure p tot in the Fluent is: That is, the total pressure is the sum of static pressure and dynamic pressure. It can be seen from Fig. 13g that the total pressure is negative, about 3500 Pa, which is consistent with the static pressure of the negative pressure cavity. It can be seen from Fig. 13h that the impeller negative pressure cavity is divided into three areas: the impeller inlet area with a smaller total pressure, the impeller outlet area with enormous dynamic pressure, and the total pressure are between the transition zone between the inlet and the outlet. It can be seen that for fluids with low flow rates, static pressure plays a significant role in total pressure; for fluids with relatively high flow rates, dynamic pressure plays a significant role in total pressure.
Analysis of simulation results of impeller and negative pressure cavity. Analysis of simulation results of negative pressure cavity. The change of the negative pressure value at the inlet of the fluid domain of the negative pressure device is shown in Fig. 14. When the number of iteration steps is 1, the calculated inlet negative pressure value is 1072.03 Pa, and then the negative pressure value rises rapidly. When the number of iteration steps is about 50, the negative pressure value reaches a stable state, 3534.75 Pa.
Under the same conditions, multiple simulations were performed on the fluid domain of the negative pressure device, thereby reducing accidental errors caused by the simulation and calculating the negative pressure value of the negative pressure cavity. The simulation results are shown in Table 4: Table 5 lists the negative pressure values of the experimental results and CFD simulation results, and the simulation results are averaged at 3534.75 Pa. For the standard sealing method, the simulation results are consistent with the experimental results, and the error is only 0.23%, which proves the correctness of the simulation www.nature.com/scientificreports/ calculation results; for the thickened sealing method, the error between the simulation and the experimental results is 5.12%, which is within an acceptable range within. In summary, CFD simulation and actual experiment have a good fit, and CFD simulation can be used instead of experiment.
Analysis of the result of impeller simulation total pressure. For an impeller whose outlet is discharged toward the atmosphere, its total pressure equals the total inlet pressure. Table 6 compares the total pressure of the three types of impellers. The first method is to calculate the total pressure difference between the outlet and the inlet of the impeller; the second method is to calculate the total inlet pressure; the third method is to measure it through experiments.
The comparison shows that the total pressure values obtained by the three methods are not much different, which shows that the CFD simulation results are consistent with the experimental results. However, there is a big gap between the actual total pressure of the impeller and the total design pressure of the impeller. The first is because the surface of the impeller and the blades are relatively rough, and the friction with the air during rapid rotation causes a pressure drop. Second, the design of the impeller is based on an empirical formula. In the process, many parameters are selected based on experience. There is a particular error, which causes a large gap between the total pressure of the impeller and the actual one.

Optimized design of impeller
Selection of impeller optimization parameters. The shape of the blade, the size of the inlet of the impeller and the number of blades have an important influence on the performance of the impeller. Therefore, the blade outlet angle, the impeller inlet diameter, and the number of blades are selected as design variables. Table 7 shows the value ranges of the three design variables. Experimental design and model selection. Before optimizing the impeller parameters, it is necessary to design them and establish a three-dimensional model for CFD simulation. In order to reduce the number of tests, a specific test method is usually selected for experimental design, which can achieve the expected results   37 .
In order to select appropriate test sample points, the Design-Expert software is used for the relevant design of the BOX method. The Design-Expert software is software specially used for test design and related data analysis. When using it for data analysis, the primary process is as shown in Fig. 15 [38][39][40] . Table 8 shows the sample points generated by the Design-Expert software, with three influencing factors and one dependent variable.
Optimal design of impeller parameters. Establish the Kriging model of the impeller. Since there is no apparent functional relationship between the negative pressure value p of the impeller evaluation index and the design variables of the impeller, the Kriging model is used to establish the mathematical model of the impeller parameters. The kriging model is a mathematical model with high calculation accuracy, small calculation amount, and short calculation time, which is very suitable for application in optimization. Simpson 41 and others have proved that the Kriging model has the characteristics of unbiased estimation and nonlinear approximation of sample points. Therefore, the Kriging model is an ideal model for the nonlinear model of the impeller 42 .
The Kriging model is a parameter difference method composed of polynomials and random process parts: Among them, y(x) is an unknown function of interest to the experimenter, T is a function of x , β = [β 1 , β 2 , ..., β n ] T is a correlation coefficient, Z(x) is a random error function with a mean value of 0, and a variance of not 0, that is: And its covariance is: Among them, R x i , x j is a n × n correlation matrix that is symmetrical along the diagonal, and R x i , x j is the correlation function between any two observation points x i and x j on the diagonal elements. The related function R x i , x j is:   where y = y 1 , y 2 , ..., y n T is the response value of f (x) . When f (x) is a constant, p is a column vector containing components n 0 , and: By maximizing the following expression, the maximum likelihood estimate θ k in the equation can be obtained: Use the sample data in Table 8 to establish a kriging surrogate model between impeller design parameters and negative pressure values. The negative pressure value corresponding to the original design parameters is predicted using the established kriging model and compared with the original data. The comparison results are shown in Table 9. The predicted value of negative pressure value of the kriging model is 3498.55 Pa, and the error with the simulation result is only 1.02%, which proves the correctness of the impeller kriging model.
Genetic algorithm optimization. Genetic algorithm (GA) is an optimization algorithm based on the evolutionary theory of "natural selection by nature, survival of the fittest." It mainly involves operations such as selection, crossover, and mutation. The higher the fitness of the individual, the greater the probability of reproducing offspring. On the contrary, if the individual's fitness is low, the probability of reproducing offspring is low, or even weeding out offspring without reproducing. Therefore, the optimal generation can be obtained through such continuous genetic optimization 43 .
According to the CFD simulation results, the negative pressure value p of the negative pressure cavity is selected as the evaluation index of the impeller parameters. The larger the negative pressure value, the better, so the objective function of impeller optimization can be expressed as: Using GA to optimize the design of the established optimization model, the value of the objective function and the optimization process of the three design variables are shown in Fig. 16.
It can be seen from Fig. 16 that the convergence value of the objective function is − 4491.19, the convergence value of the blade outlet angle is 54.55°, the convergence value of the impeller inlet diameter is 67.43 mm, and the convergence value of the number of blades Z is 18. The optimized result is the same as the original design. The comparison is shown in the Table 10. Compared with the original design, the negative pressure value has been increased from 3534.75 Pa to 4491.19 Pa, increasing 27.06%. Therefore, the optimized design results of the impeller parameters are: the blade outlet angle is 54.55°, the impeller inlet diameter is 67.43 mm, and the number of blades Z is 18. Under this set of parameters, the corresponding negative pressure value is 4491.19 Pa.

Conclusion
This paper takes the negative pressure adsorption wall-climbing robot as the research object. Through CFD simulation and prototype experiment, the influence of wall material, adsorption position, shape of the negative pressure cavity and sealing method on the adsorption performance of the wall-climbing robot is studied, and it is based on Kriging. The model and genetic algorithm optimize the parameters of the impeller in the adsorption device, and the following conclusions are obtained: 1. Without considering the deformation of the wall surface, the wall material has little effect on the adsorption performance. 2. By comparing the negative pressure value when the robot is adsorbed in three different poses, it can be considered that the posture of the wall-climbing robot does not affect the negative pressure value.   Figure 16. Optimization of objective function and design variables.